arXiv: 1505.04192v2 [hep-ph] 15 Jan 2016 


UMD-PP-014-018 


A Hidden Dark Matter Sector, Dark Radiation, and the CMB 

Zackaria Chacko, 1 Yanou Cui, 1,2 Sungwoo Hong, 1 and Takemichi Okui 3 

1 Maryland Center for Fundamental Physics, Department of Physics, 

University of Maryland, College Park, MD 20742, USA 
2 Perimeter Institute for Theoretical Physics, Waterloo, Ontario N2L 2Y5, Canada 
3 Department of Physics, Florida State University, Tallahassee, FL 32306, USA 

We consider theories where dark matter is composed of a thermal relic of weak scale mass, whose 
couplings to the Standard Model (SM) are however too small to give rise to the observed abundance. 
Instead, the abundance is set by annihilation to light hidden sector states that carry no charges under 
the SM gauge interactions. In such a scenario the constraints from direct and indirect detection, 
and from collider searches for dark matter, can easily be satisfied. The masses of such light hidden 
states can be protected by symmetry if they are Nambu-Goldstone bosons, fermions, or gauge 
bosons. These states can then contribute to the cosmic energy density as dark radiation, leading 
to observable signals in the cosmic microwave background (CMB). Furthermore, depending on 
whether or not the light hidden sector states self-interact, the fraction of the total energy density 
that free-streams is either decreased or increased, leading to characteristic effects on both the scalar 
and tensor components of the CMB anisotropy that allows these two cases to be distinguished. The 
magnitude of these signals depends on the number of light degrees of freedom in the hidden sector, 
and on the temperature at which it kinetically decouples from the SM. We consider a simple model 
that realizes this scenario, based on a framework in which the SM and hidden sector are initially in 
thermal equilibrium through the Higgs portal, and show that the resulting signals are compatible 
with recent Planck results, while large enough to be detected in upcoming experiments such as 
CMBPol and CMB Stage-IV. Invisible decays of the Higgs into hidden sector states at colliders can 
offer a complementary probe of this model. 

I. INTRODUCTION 

In the last two decades, with the advent of precision 
cosmology, it has become clear that some form of non- 
luminous dark matter (DM) contributes more than 20% 
of the total energy density of the universe Q], Although 
it is known that the particles of which DM is composed 
lie outside the SM of particle physics, their precise nature 
remains to be understood. 

In the absence of a detailed understanding about the 
properties of dark matter, many different candidates have 
been put forward. A large class of well-motivated theories 
are based on the ‘Weakly Interacting Massive Particle’ 

(WIMP) paradigm. In its simplest incarnation, this sce¬ 
nario involves a particle of weak scale mass, the WIMP, 
that has interactions of weak scale strength with the SM 
fields. This class of theories possesses the very attractive 
feature that the WIMPs that survive after their anni¬ 
hilation into SM particles freezes out naturally tend to 
have a relic abundance that is in good agreement with 
observations M- 

In this conventional scenario, the WIMP must have 
interactions of weak scale strength with the SM fields. 

Several different types of DM experiments are searching 
for evidence of such interactions. These include direct 
detection experiments that are looking for the recoils of 
nuclei after being impacted by a DM particle, indirect 
detection experiments that seek to observe the products 
of DM annihilation, and collider experiments such as the 
the Large Hadron Collider (LHC) that seek to produce 


DM. Till date, there has been no compelling evidence for 
the existence of such interactions, and the experimental 
limits now exclude a significant part of the preferred pa¬ 
rameter space for many WIMP DM candidates HI]. 

With the simplest realizations of the WIMP paradigm 
beginning to come under strain, several ideas have been 
put forward to explain the absence of a signal in these 
experiments. Among the hypotheses that have been 
advanced are that the DM candidate scatters inelasti- 
cally HUD], is leptophilic CMS], or interacts preferen¬ 
tially with heavier quark flavors CMH]- An alternative 
proposal that has attracted interest [lSHSS] is the idea 
that, while DM is indeed composed of WIMPs, their cou¬ 
plings to the SM fields are suppressed, and too small to 
yield the observed abundance. Instead, the DM candi¬ 
date possesses interactions of weak scale strength with a 
new hidden sector that carries no charge under the SM 
gauge interactions, and its relic abundance is set by its 
annihilation into these states. Such a scenario can nat¬ 
urally account for the observed abundance of DM, while 
explaining the absence of any signal in experiments. 

It is not difficult to envisage scenarios where the DM 
candidate naturally has weak scale mass and interactions 
of weak scale strength with a hidden sector. For exam¬ 
ple, in supersymmetric theories, both the weak scale and 
the scales in the hidden sector could be set by the scale 
of supersymmetry breaking. Similarly, in extra dimen¬ 
sional Randall-Sundrum constructions, both the Higgs 
and the hidden sector states could be localized to the 
infrared brane. In such a scenario the scales in the hid- 
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den sector would again naturally be of order the weak 
scale. Therefore, provided the SM and hidden sectors 
are in thermal equilibrium at or above the weak scale, 
so that their temperatures at freeze out are not very dif¬ 
ferent, this framework can naturally explain the observed 
abundance of DM while remaining consistent with all ex¬ 
perimental constraints. 

The existence of a hidden sector into which DM anni¬ 
hilates can potentially be tested by experiments. The na¬ 
ture of the signals depends on the masses of the particles 
in the hidden sector, and on their couplings to SM states. 
If all the particles in the hidden sector have masses above 
an eV, and the temperature of this sector is comparable 
that of the SM, we expect these states will decay or an¬ 
nihilate into SM particles before the CMB epoch. This 
is because the lightest state in the hidden sector, be¬ 
ing massive, would otherwise contribute to the energy 
density in DM, violating the overclosure bounds if it is 
heavier than a keV, and coming into conflict with the cos¬ 
mological constraints on a warm sub-component of DM 
if it is lighter than a keV. Such a scenario therefore im¬ 
plies the existence of couplings between the hidden sector 
states and the SM that can potentially be tested in ex¬ 
periments, as in the scenarios of exciting DM [19: 23] . 
secluded DM [20] and boosted DM [23] [25]. If instead, 
some or all of the states in the hidden sector have masses 
below an eV, they would be expected to constitute a sig¬ 
nificant component of the energy density of the universe 
both before and during the epoch of matter-radiation 
equality, potentially leading to observable signals in the 
CMB (2TJ 1221 I2GH31| . The simplest possibility is that 
these states, if present, are massless, and constitute dark 
radiation (DR) at present times, thereby obviating the 
need for any other mass scales in the theory. It is this 
scenario, and the associated signals, that we will focus 
on in this paper Q 

The presence of these new light particles implies the 
existence of additional structure in the theory, if the 
scenario is to be natural. There are three known sym¬ 
metries that can prevent masses from being generated 
for a massless particle: a shift symmetry for a spin-0 
Nambu-Goldstone boson, a chiral symmetry for a spin- 
1/2 fermion, and a gauge symmetry for a spin-1 vector 
boson. DR candidates protected by these symmetries 
have been considered, for example, in [3H1 [37]. These 
symmetries may appear individually or in combination; 
for example, the spectrum of light states may consist of 
a single Nambu-Goldstone boson, but it may also consist 
of spin-1/2 fermions charged under a U(l) gauge group 
with its associated massless spin-1 boson [38]. This latter 


1 If the DM and DR are tightly coupled, oscillations of the DM- 
DR fluid can also give rise to signals in the matter power spec¬ 
trum |321 188 1 . (see also |34l I35 | h However, the large strength 
of the interaction required to obtain an observable effect would 
overly deplete the abundance of DM, and is therefore disfavored 
in the framework of thermal WIMPs. 


example illustrates that the constituents of the DR need 
not be free, but may have interactions amongst them¬ 
selves without violating the symmetries that keep them 
light. In general, therefore, we see that the DR can take 
two distinct forms: 

• Free DR, which free streams during the era of 
acoustic oscillations, and is characterized by a mean 
free path H~ x , where H is the Hubble constant. 

• Scattering DR, which scatters during the era of 
acoustic oscillations, and is characterized by a mean 
free path <C H~ l . 

Since the presence of DR is a robust prediction of this 
scenario, it is important to understand whether it can be 
detected, and whether we can distinguish between the 
two different cases of free DR and scattering DR. It 
is these questions that we shall be primarily concerned 
with in this paper. We find that, provided the hidden 
DM sector was in thermal equilibrium with the SM at 
temperatures at or above the weak scale, the contribu¬ 
tion of the DR to the energy density during the CMB 
epoch is in general large enough to be detected in future 
experiments, such as CMBPol f32] {vN Bil = 0.044), and 
eventually CMB Stage-IV |30[ =0.02). 

We also find that it is, in general, possible to dis¬ 
tinguish between scenarios with free streaming DR and 
scattering DR. Studies of the scalar SM3] and ten¬ 
sor [331 [45] metric perturbations have established that 
the details of the CMB spectrum depend on the fraction 
of the energy density in radiation that is free streaming. 
This ratio impacts not just the amplitudes of the modes, 
but also the locations of the peaks in the CMB spectrum. 
In scenarios where neutrinos scatter off new light states 
during the period immediately prior to the CMB epoch, 
as in models of late time neutrino masses (3511331 , and in 
the neutrinoless universe scenario [49], this ratio differs 
significantly from the SM prediction. Consequently, it 
has been possible to establish that this class of theories, 
which was already disfavored by the WMAP data [13- 
[54], is now excluded by Planck, unless the new neutrino 
interactions come into equilibrium only very shortly prior 
to matter-radiation equality [55,155]. 

Similar considerations apply to the class of theories we 
are considering. In the presence of a new dark component 
of radiation, the free streaming fraction is altered, with 
the sign of the correction dependent on whether the DR 
scatters or free streams. Consequently, the amplitudes of 
the scalar and tensor modes receive corrections, with the 
sign of the effect dependent on whether the DR is free or 
self-interacting. In addition, the locations of the CMB 
peaks are shifted, with the sign of the shift again depen¬ 
dent on whether or not the DR carries self-interactions. 
We find that these effects may be large enough to allow 
upcoming experiments to distinguish between free DR 
and scattering DR. Therefore the CMB offers a window 
into the dynamics of the hidden sector that DM annihi¬ 
lates into. 
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The outline of this paper is as follows. In the next 
section we discuss the CMB signals associated with DR, 
and explain how scenarios with free streaming DR and 
scattering DR can be distinguished. In section III we 
show how the expression for the relic abundance of DM 
may be generalized to the case when the temperature of 
the dark sector differs from that of the SM. In section 
IV we consider a simple model based on the Higgs por¬ 
tal that realizes this scenario, and show that the signals 
can be large enough to be detected by upcoming CMB 
experiments, while remaining compatible with the recent 
Planck results [T]. We also show that invisible decays of 
the Higgs into hidden sector states at colliders can offer 
a complementary probe of this scenario. Our conclusions 
are in section V. 


II. THE CMB SIGNALS OF DARK RADIATION 

The CMB spectrum is affected by the presence of extra 
relativistic degrees of freedom during the era between 
matter-radiation equality and photon decoupling. It is 
customary in cosmology to quantify the contribution to 
the energy density from such additional radiation in unit 
of the energy density of a single relativistic SM neutrino 
species p w , 


AN ef f = 


_ PdR 


( 1 ) 


where p OR is the energy density of DR, and all SM neutri¬ 
nos are treated as being relativistic at the temperatures 
in question. 

For any specific hidden sector model, we can calculate 
p DR , and hence A N e g. The first step is to determine the 
temperature of the dark sector, X, that corresponds to a 
given SM temperature X at the same cosmic time t. To 
do this, note that the comoving entropies of the SM and 
of the dark sector are separately conserved after the two 
sectors thermally decouple from each other. Then, after 
thermal decoupling when X = X = Xkd, but before the 
neutrino decoupling, taking the ratio of the two entropy 
conservation relations leads to the relation 


9* T g*kd 

9 * 9* kd 


Here g* and g*kd are the number of degrees of freedom 
in the SM at temperatures X and X^d respectively, with 
the usual 7/8 factors for the fermions. The corresponding 
parameters in the dark sector are labelled by g* and g*kd- 
Applying this relation just above the neutrino decoupling 
temperature X ~ 0(10) MeV, g* = 10.75, and noting 
that the contribution of a single neutrino species to the 
energy density is given by p lv = ^fo X > we have 


AIV eff 


g*T 4 _ 4 /g* ff*kd Y /3 

7 T 4 - 7 g* g*kd / 


(3) 


Note that the above A N e g computed for the time just 
before neutrino decoupling is the same as the AN e g at 
the later CMB time, as X and T v redshifts the same way 
till then. 

As outlined in section [T[ in general p DR can consist of 
two qualitatively very different types of radiation: free- 
streaming radiation with mean free path H~ 4 , and 
scattering radiation with mean free path -C H ~V We 
can parametrize each of these components of radiation 
in complete analogy with the definition ([!]), 


D 

AA* ee = 

PlU 


AJV e ff att = 


„scatt 
P DR 


Piv 


(4) 


so that the total extra radiation AN e g = AJV e j ee + 
AlV|g att . 

In this class of models, for a given thermal decoupling 
temperature Xkd, there is a robust lower bound on AIV e ff. 
To understand this, note that the lowest possible value 
of g* in Eq. (J3| is 1, corresponding to the case when the 
dark radiation consists of just a single real scalar. Then, 
if thermal decoupling between the hidden sector and the 
SM occurs at temperatures well below the mass of the 
DM particle, we can have g*kd = g* = 1. In this limit we 
obtain a lower bound on AN e g, 


A / \ 4 /3 

AA e t=z • (5) 

' \9*kdJ 

From Eq. ([5|, assuming all the SM degrees of freedom 
are in the bath at decoupling, we have g*kd = 106.75, 
which leads to a lower bound on the effective number of 
neutrinos, AN^ n > 0.027. This result applies to arbi¬ 
trarily high Xkd provided there are no new states in the 
SM sector up to that scale. In Fig.[l]we have plotted this 
lower bound as a function of the decoupling temperature 
Xkd- Quite intriguingly, the ultimate experimental sensi¬ 
tivity at CMB-Stage-IV is <7n bH = 0.02 [30], which would 
allow it to probe the DM scenario we outline here. This 
projected experimental sensitivity is based on a one pa¬ 
rameter extension of the standard six parameter ACDM 
model that accommodates varying Neg. Although this 
projection assumes that the primordial density fluctua¬ 
tions have an exactly power law spectrum, the recent 
Planck results have established that at present this is an 
excellent fit to data [T], and so this forms a reasonable 
basis for estimating the sensitivity. 

If the dark sector consists of just a real scalar, we ex¬ 
pect that it is a Goldstone boson, so that its mass is pro¬ 
tected against radiative corrections from the weak scale. 
In this scenario, the interactions of the states that consti¬ 
tute the DR are momentum suppressed, and so the DR 
free streams. In scenarios where the DR possesses self 
interactions large enough to prevent free streaming, the 
requirement of naturalness up to the weak scale implies 
that it must be composed of more than just a single real 
scalar, or else the radiative corrections to the scalar mass 
from the self interactions would tend to make its mass 
much greater than an eV. Therefore, in scenarios where 
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FIG. 1. A N e ff as a function of the temperature at which 
the SM and dark sector thermally decouple. Also shown are 
the 2015 Planck results: the central value (Green dashed line) 
and the 2<r constraint (Orange dashed line). 


the DR scatters, we expect that there will be additional 
light states in the hidden sector, and so A N e g is expected 
to be larger than its minimum value, A N£g n . In Fig. |TJ 
we have plotted AN e g as a function of Tkd for the case 
when the DR consists of a pair of massless Weyl fermions 
with vector-like charges under a U(l) gauge group, and 
the associated massless gauge boson. We see that even 
for high Tkd, we predict AlVeft > 0.15, which is large 
enough to be observed at CMBPol. 

In the rest of this section we discuss how the exper¬ 
imental limits on A N e g are obtained, and how we can 
distinguish between the two cases of scattering DR and 
free streaming DR. 


A. The Determination of AA e fr 

At present, limits on AN e g are obtained by considering 
how the presence of additional energy density in radia¬ 
tion would affect the quality of the fit in the six parameter 
ACDM model. Of the six parameters, two are particu¬ 
larly sensitive to A N e g. These are the total energy den¬ 
sity in matter, p m , which is the sum of the energy densi¬ 
ties in baryons and DM, and in the cosmological constant, 
Pa- The presence of additional energy density in radia¬ 
tion would tend to delay the onset of matter-radiation 
equality. Since the amplitude of a Fourier mode is very 
sensitive to the fraction of energy density in matter as 
it crosses the horizon, this ratio is highly constrained by 
the CMB data. Therefore, for AlV e ff > 0, the best fit 
is obtained by increasing p m in the appropriate propor¬ 
tion to ensure that the redshift at the onset of matter- 
radiation equality is unaffected. Since the energy density 
in baryons, pb , is very tightly constrained by measure¬ 
ments of the relative heights of the even and odd CMB 
peaks and cannot be altered, the change in p m is accom¬ 
plished by an increase in the energy density in DM. 


The additional energy density in DR, and in matter, 
then implies an increase in the Hubble constant during 
the CMB epoch. This will in turn affect the size of the 
sound horizon, leading to a change in the locations of 
the CMB peaks. This observable is, once again, highly 
constrained by the data. However, this effect can be off¬ 
set by changing pa so as to alter the distance to the last 
scattering surface, thereby keeping the angular locations 
of the peaks intact. Nevertheless, as we now explain, the 
change in the Hubble constant during the era of acoustic 
oscillations leads to other effects in the CMB spectrum 
that can no longer be compensated for once p m and pa 
are fixed. 

Prior to recombination, the photons interacted 
strongly with the baryons. Although the photon mean 
free path during this era was relatively short, the pho¬ 
tons were nevertheless able to diffuse outward, with a 
characteristic diffusion length r^. As a consequence of 
this diffusion, inhomogeneities and anisotropies at scale 
smaller than the are suppressed. This damps the 
peak amplitudes at higher £ relative to the first peak 
at £ ~ 220, which corresponds to modes that entered the 
horizon near recombination. This effect is known as Silk 
damping, or diffusion damping. A change in the Hub¬ 
ble rate affects the time available for diffusion, leading 
to observable effects on the CMB spectrum. In particu¬ 
lar, the height of the first CMB acoustic peak relative to 
the latter peaks is altered. Therefore, this effect can be 
used to place limits on the Hubble constant during the 
epoch of acoustic oscillations, and therefore on A N e g. 
The presence of additional energy density in radiation 
also leads to changes in the CMB spectrum associated 
with the early Integrated Sachs-Wolfe (ISW) effect, but 
these are less significant than the effects arising from Silk 
damping El- 

In principle, an increase in the fraction of baryons in 
helium, Yne, while pb is held fixed, would reduce the num¬ 
ber of free electrons available for scattering, and could 
also account for a change in the scale of Silk damping. 
However, the helium fraction in the SM can be calculated 
sufficiently precisely from Big Bang nucleosynthesis so as 
to exclude this as the explanation for any observed dis¬ 
crepancy. For a good discussion of these issues with more 
details, see mm 

B. Distinguishing between Free and Scattering DR 
via Scalar Metric Perturbations 

Several authors have considered the effects of the SM 
neutrinos on the scalar component of the CMB spec¬ 
trum mm\ . These results can easily be generalized 
to the case when there is additional energy density in 
radiation, and can be used to distinguish between free 
streaming DR and scattering DR. 

In the conformal Newtonian gauge the Robertson- 
Walker metric with scalar perturbations takes the form, 

ds 2 = a 2 (r) (- (1 + 2$) dr 2 + (1 - 2tf) dr 2 ) (6) 
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Here r represents conformal time, while a is the cosmo¬ 
logical scale factor, T and $ represent the scalar metric 
perturbations. In the absence of any free streaming par¬ 
ticle species, we have 4' = 4>. When, however, a free 
streaming species is present, the energy momentum ten¬ 
sor becomes anisotropic. This leads to a difference be¬ 
tween and 4> that is proportional to /„, the total energy 
density in free streaming radiation expressed as a fraction 
of the total energy density in radiation. 


u = 


Pall free rad 
Pall rad 


3pi 


, free 
1 rDR 


3pi 


Pi 


■ 


pscatt 


(7) 


In the limit that and Por U are small compared to 
p a ii ra( j, the total energy density in radiation, the devia¬ 
tion from the standard cosmology is given by 


fu ~ f, 

f, 


SM 

SM 


( 1 -/,|sm) ANlr-UlsuAN^ 


3 

= (0.59AIVgg ee - 0.4lAA|g att ) . 


( 8 ) 


Now, the solution of the coupled system of equations 
for matter, radiation and gravity reveals that the pres¬ 
ence of a free streaming component in radiation is as¬ 
sociated with a change in the amplitudes of the CMB 
modes at large £. The magnitude of this effect was first 
determined numerically in [41] . Subsequently, analytic 
expressions were obtained in [42l 03] . The result is given 
by, 


SCe 

Ci 15 U 


(9) 


Then, using Eq. ([8]), we can obtain an expression for 
the fractional change in Ci with respect to the standard 
cosmology, 


A C t 
C t 


SC, 6C t 


Ct 


Ci 


SM 


15 v" I sm) 

—0.072 (0.59AA e f ^ ee - 0.41AA e s £ att ) . 


( 10 ) 


We see that the result is independent of l, and that the 
sign of this effect depends on whether the DR is scatter¬ 
ing or free streaming. 

In addition to the corrections to the amplitude, there 
is a shift in the angular locations of the high i CMB 
peaks by an equal amount [43] . This signal is particularly 
important because, in contrast to other effects of DR such 
as Silk damping, it is difficult to mimic by altering other 
parameters such as the helium fraction. The magnitude 
of this shift is again proportional to the free streaming 
fraction /„, 


6£ 



( 11 ) 


Here t\ ~ 300 represents the average angular spacing be¬ 
tween the CMB peaks at large £. Again, in the limit that 
DR contributes only a small fraction of the total energy 
in radiation, AAT e ff < 1, Eq. Q leads to an expression for 
the change in Si with respect to the standard cosmology, 


A£ = 5l-5i\ SM 


= -57(/„-/„| sm )- 

~ -7.8 (0.59A 


£a_ 

300 

- 0.41AA e s £ att ) 


£A_ 

300 ‘ 


( 12 ) 


Once again we see that the sign of the effect depends on 
whether the DR is scattering or free streaming. 

We can obtain a very rough estimate of the sensitiv¬ 
ity of upcoming CMB experiments to the effects of At 
by considering how well N e g- can be determined when the 
helium fraction Yho is allowed to float freely. In this limit, 
the effects of A N e g on Silk damping can be compensated 
for by changes in F He . Under these circumstances, the 
shifts in the locations of the CMB peaks play an im¬ 
portant role in the determination of N e ff, and we can 
interpret the results as a rough guide to the sensitivity of 
these experiments to A N e g arising from its effect on At, 
and not its effect on Silk damping. The projected sen¬ 
sitivity of CMBPol to N e ff when Ifje is allowed to float 
is A N e g = 0.09 [43]. We therefore expect that provided 
A N e g > 0.10, upcoming experiments will have some sen¬ 
sitivity to whether DR is free streaming or scattering, 
allowing the possibility of distinguishing between these 
two scenarios. 


C. Distinguishing between Free and Scattering DR 
via Tensor Metric Perturbations 

The presence of a free streaming component of radia¬ 
tion also affects the tensor component of the CMB spec¬ 
trum. Detailed studies of the effects of the SM neutri¬ 
nos on the tensor modes (the B- and Ill-modes) of the 
CMB were performed in [10105], which found an 0(10)% 
damping of the correlation functions of the tensor modes 
at long wavelengths, rising to an 0(35)% damping at 
short wavelengths. Analytic results for the damping were 
subsequently obtained in ISSUED]. The corrections to the 
spectrum that arise from the presence of a free streaming 
DR component were considered in [61] . 

We now show that the results of [44] can be generalized 
in a very simple way to arbitrary N e g, provided AN e g < 
1. The crucial observation is that, in the analysis of 144] . 
the effects of the SM neutrinos arise entirely from their 
contribution to /„, the free-streaming fraction of the total 
energy density, 

j _ Pall free rad _ ^Plu T Pdr (13) 

Ptotal Ptotal 

Therefore, by understanding how the result depends on 
/„, we can immediately determine how the correlation 
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functions of the tensor modes depend on AN^ e and 
AfV|g att . During the radiation dominated era, to a very 
good approximation, f v = /„. However, as matter- 
radiation equality approaches, the contribution of matter 
to the total energy density can no longer be neglected, 
and /„ and f v are distinct. 

The Robertson-Walker metric with tensor perturba¬ 
tions takes the form, 

ds 2 = a 2 (r) (—dr 2 + [<5^- + /i,y(x, r)] dx l dx J ) (14) 

with ha = dihij = 0. We define a new coordinate 
u = hr, where k is the co-moving wave number. Then, 
as shown in [44] , the amplitude of tensor perturba¬ 
tions with co-moving wave number k can be written as 
hij(u) = hij( 0) x( u )> where the function \{u) remains to 
be determined. It satisfies the integro-differential equa¬ 
tion 

JY{u,d u )x{.u) = f v [ l{u, U)x'(U)dU, (15) 
Jo 

where T(u,d u ) is a second-order, linear differential op¬ 
erator. Its precise form, along with that of the kernel 
function I(u,U), may be found in [33] • The initial con¬ 
ditions on x are given by x(0) = 1 and x'(0) = 0. 


In general, the integro-differential equation Eq. (15) is 


not simple to solve. However, approximate analytic so¬ 
lutions that apply in certain limits have been obtained 
for the important case of the SM with 3 free streaming 
neutrinos, N e s = 3.046. As we now explain, these solu¬ 
tions can be generalized in a simple way to the case of 
arbitrary /„, provided AlV e ff < 1. 

For short wavelengths that enter the horizon well be¬ 
fore matter radiation equality, that is, u 1, the solution 
to Eq. (15) approaches a homogeneous solution [44] . 


X{u) ->• A 


sin(u + <5) 


where the parameters A and S contain the dependence 
on fv = fv In the limit /„ = 0, A and <5 take values 
A 0 = 1, <5o = 0- We denote this solution by xo(w). For 
the case of the three free streaming neutrino species of the 
SM, with fv = fv M = 0.41, a numerical study [H] leads 
to the values Asm = 0.80, <5sm ~ 0. The fact that the 
value of A changes by only about 20% for the change from 
iV ef j = 0 to N e f[ = 3.046 indicates that for A N e g < 1 the 
term proportional to f v can be treated as a perturbation. 
Accordingly, we may obtain an approximate solution for 
general /„ by replacing x( u ) by Xsm(w) on the right hand 
side of Eq. (15). Here Xsm(w) is the solution of Eq. (15) 
for fv = Having made this approximation, Eq. (15) 
reduces to 

d u ) X (u) = fv f l(u, U) x'sm (U) d U . (17) 

Jo 

Noting that the right hand side of this equation is pro¬ 
portional to J-(u, d u ) xsm(u), we have 


_u_ 

JSM 


XSM (u) 1=0. 


(18) 


Recalling that F(u,d u )x o(w) = 0, it follows that this 
equation admits a solution of the form, 


fv 

x(«) - JsmXsm(u) oc xo (u) 


(19) 


Using the initial conditions to fix the constant of propor¬ 
tionality, we obtain an analytic solution for x( u ) valid 
for general f v that is applicable in the short-wavelength 
limit, 

xM^fl + i^SM-l)) — ■ (20) 


J l 


As discussed in [ 33 ] , the ratio of the CMB tensor corre¬ 
lation functions for general f v relative to that for f 8M is 
given by 


R = 


X'(u) 


Xsm(' 1 


1 + 


FSM 


(Asm — 1) 


( 21 ) 


A 2 

^SM 


Since the right hand side of this equation is indepen¬ 
dent of u, it follows that the fractional change in the 
correlation functions arising from the presence of DR is 
independent of wave number for modes deep inside the 
horizon. The result is plotted in Fig. 2. 

To obtain the solution of Eq. (15) for long wave¬ 


lengths , it is convenient to change variables from u to 
y = a(r)/a(r eq ), where a(r) is the standard scale fac¬ 
tor of the Robertson-Walker metric and r eq is the value 
of the conformal time coordinate r at matter-radiation 
equality. Expressed in terms of y, the integro-differential 
equation Eq. (15) takes the form 


(16) [S(y, d y ) + Q 2 ] x(y) = fv r J(y, Y) x'(U) du, (22) 

Jo 


Here Q(y,d y ) is again a second-order, linear differential 
operator, and Q is a normalized comoving wave number 
defined as Q = \/2k/k eq , where k eq is the value of k 
corresponding to the length scale that enters the horizon 
at matter-radiation equality. We have also eliminated 
fv{y) in favor of fv, which is independent of y. The initial 
conditions are now given by x(0) = 1 an( l X 7 (0) = 0. The 
precise forms of Q and the kernel function J(y , Y) may be 
found in |44j . Then, given the solutions Xo(y) and Xsm(i/) 
corresponding to the choices f v = 0 and f v = /® M (y), an 
approximate solution x(u) corresponding to a general f v 
may be obtained as 

x(y) = xo(y) + j^m(xsm( y) - xo(y)) ■ (23) 

Expressions for XsM^dec) and x'oiVdec ) at small Q were 
obtained in m, 

X S M(?/dec) = a 2Q 2 + O 4 Q 4 + 06 Q 6 + 0(Q 8 ) , 

Xo(2/dec) = b 2 Q 2 + b 4 Q 4 + b 6 Q 6 + 0(Q 8 ). 


-F(m,<9„) x(u) 


(24) 
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Tensor mode anisotropy (short wavelengths) 



FIG. 2. Relative change in tensor mode anisotropy compared 

/ / \ 2 
x («) 


to the SM case (N e g = 3.046), R = 
wavelengths. 


Xsm( u ) 


for short 


Tensor mode anisotropy (long wavelengths) 



- SM (an fs =o, an sc =o) 

-an fs =i,an sc =o 

- AN fs =0.5, AN S c=0 

. AN fs = 0.2, AN sc =0 

-AN fs =0, ANsc=1 

-AN FS =0, AN S c=0.5 

. AN fs =0, AN sc =0.2 


FIG. 3. Relative change in tensor mode anisotropy compared 


to the SM case (N e g = 3.046), R = 


X (Vdec) 


AN^ e (shorthand AN FS ) and A N°g 
for long wavelengths. 


, for varying 

XsM (?/dec) 

(shorthand ANsc), 


with 

o 2 = -0.573661, a 4 = 0.243294, a 6 = -0.0381643, 
b 2 = -0.601254, 6 4 = 0.264482 , b 6 = -0.0424186 . 


Here ?/dec is the value of y at photon decoupling, given 
by t/dec = 3.31. Note that the regime of validity of these 
expansions is limited to small Q, and they are expected 
to break down at Q ~ 1. Nonetheless these analytic 
approximations offer a simple way to parametrize the 
effects of DR on the tensor mode anisotropies for long 
wavelengths. Using these expressions in association with 
the general solution Eq. (23), we can determine the sup¬ 
pression of the tensor correlation function for arbitrary 
f v relative to that for the SM, 


X’iVdec) 2 

XsM (pdec) 


(25) 


The difference between the scenario we are considering 
and the conventional relic abundance calculation for a 
thermal WIMP lies in the fact that in Eq. (26), the aver¬ 
aged annihilation cross section (cr v) and the DM equilib¬ 
rium number density n® q now depend on the dark sector 

temperature T rather than the SM temperature T, while 
the Hubble constant H depends on both T and T. In 
our analysis we will assume that at temperatures in the 
neighborhood of the freeze-out, the number of relativis¬ 
tic degrees of freedom in the SM, and in the dark sector, 
are not changing. Then, since the entropy densities of 
the two sectors are separately conserved, it follows that 
during freeze-out the ratio of the dark sector tempera¬ 
ture T to the SM temperature T is a constant labelled 
by r = T/T. We can then express all the terms in the 
Boltzmann equation as functions of the SM temperature 
T. For the equilibrium number density in the nonrela- 
tivistic limit we have, 


The results are shown in Fig. [3j where we can see that 
scattering and free streaming DR contribute with oppo¬ 
site signs. It follows that while scattering DR tends to 
reduce the damping effect associated with the SM neu¬ 
trinos, free streaming DR enhances this effect. This dis¬ 
similar behaviour is analogous to what was observed for 
the corrections to the amplitude and phase of the scalar 
modes. 


III. DETERMINATION OF THE RELIC 
ABUNDANCE 

In this section we explain how the well-known for¬ 
malism for determining the relic abundance of DM gen¬ 
eralizes to the case when the DM sector and the SM 
are at two different temperatures. The evolution of the 
DM density in an expanding universe is governed by the 
Boltzmann equation, 

^ + mnx = _ {av) („2 _ ( n e q)2) (26 ) 


n e x q (T)=g x 

= 9x 

For the Hubble constant, 


m x T\ 

uu);' 

( rmJT \ 2 
1 a 1 e rT . 


^ 2t r 


H = 


87 tG 


8n 


~P = 


3Mp 


(pSM + Pdark) 


3 

- 877 2 (—g*{T) T 4 + — i/*(T) f 4 
3M|; \30 y v ' 30 y v ’ 

47T 3 4 

45M| 5 * eff ’ 


where 


(27) 


(28) 


g* e ff(T, T) = g*(T) +g*(T) I - I = g*{T) +p*(T)r 4 . 


( 29 ) 











































In the nonrelativistic regime (av) is well approximated 
by 

(av) = a + b(v 2 ) (30) 


where the constant term plays the leading role if the anni¬ 
hilation can proceed through the s-wave. If, however, the 
s-wave contribution is suppressed so that annihilation oc¬ 
curs primarily through the p-wave, the term proportional 
to v 2 dominates. Performing the thermal average, 


T 


a + b(v 2 ) = a + 6 rb -. 

m v 


(31) 


We are now in a position to solve the Boltzmann equa¬ 
tion by the standard procedure. We introduce the vari¬ 
able Y = n x /ssM , where ssm is the entropy density car¬ 
ried by the SM degrees of freedom. We also introduce 
the variable x = m x /T . In the radiation dominated era, 
when g * and g* are not changing, we have 

£ = Hz (32) 

at 

Eliminating the number density n x in favor of Y, and t 
in favor of x, we obtain for the Boltzmann equation, 


dF 

dx 


(<™> „ (V 2 V 2 \ 
SSM (Y - l eq ) 


(33) 


Hx eq/ 

Rewriting this in terms of A = Y — r eq , we have 

A' = -K q — f(x) A (2F eq + A) (34) 

where f(x) is given by 

6 rb/x 


(ctu)ssm 
!lx) = Ijx 


45 yj 


9* a ' 

m x M p — 


(35) 


This Boltzmann equation cannot be solved exactly. How¬ 
ever, it is possible to obtain analytic solutions of this 
equation that are valid at very early times, and very late 
times. Then, an approximate solution that is valid for all 
times may be obtained by stitching together these lim¬ 
iting cases. We define the freeze-out temperature Tf in 
terms of the implicit relation, 


A (x f ) = cY eq (x f ) , (36) 


where c is an order one number. The final result depends 
only logarithmically on the value of c, which is chosen to 
be y/2 — 1 in the case of s-wave annihilation and \/3 — 1 
in the case of p-wave annihilation, to give the best fit to 
numerical results [B5]- At early times i<Ci/, A'< T/ q , 
so that the Boltzmann equation reduces to 


A = — 


Y' 

eq 


f(x) (2F eq + A) 


(37) 


This equation, in combination with Eq. (36), may be used 
to determine the freeze-out temperature, 


Xf = r log 


/ c(c + 2) 
y 4-7T 3 



9x 

\/ 


m x Mp 


r 2 (a + 6 rb/xf) 


1/2 


( 38 ) 


At late times, x Xf, A « F l^ q and A' F eq so 
that we have, 


Y^Y 1 = -/(*) 

Integrating this equation from Xf to 00 , we obtain 


Y = 

■ L OO 


45 yj 


Mp TO, 


a + 3rb/xf 
x f 


y -1 


(39) 


(40) 


where we have used the fact that A (xf) A ( 00 ). 
Combining this result with the expression for Xf from 
Eq. (38), we can obtain the present day energy density 
in dark matter, f l x = to^So^ooPc 1 - Here so, p c are the 
present-day SM entropy density and critical density, re¬ 
spectively. 

From this discussion we see that the difference between 
our framework, which allows for T ^ T at freeze-out, and 
the conventional scenario with T = T, primarily trans¬ 
lates into an extra factor of r = T/T in the expression for 
the relic abundance of DM. There is an additional effect 
arising from the r dependence of the argument of the log¬ 
arithm in Eq. (38), but this is small. Put in another way, 


the values of a and b in the expression for (av) in Eq. (30) 
that correspond to the observed Hdm are smaller by a 
factor of T/T than in the case of the standard thermal 
cross section ~ 3 x 10~ 26 cm 3 /s. One may obtain similar 
results to those given in Eqs. (38) and (40) by performing 


a simple estimate based on equating P x = n^(av) to the 
Hubble constant H at freeze-out, while keeping track of 
the distinction between T and T. 

In this analysis, we have assumed that the number of 
degrees of freedom in the SM, and in the hidden sector, 
do not change in the neighborhood of the freeze-out tem¬ 
perature. We have further assumed that the annihilation 
cross section is well approximated by an expansion of the 
form given in Eq. (l3l|). Using the methods of [63, !B31 . 


these assumptions can be relaxed, and the result gener¬ 
alized. We leave this for future work. 


IV. A SIMPLE BENCHMARK MODEL 
A. The Model 

In this section we consider in detail a simple model 
that illustrates the scenario we are considering. We con¬ 
sider an unbroken U(l) gauge theory with massless spin- 
1/2 fermions charged under it. Let A^ be the massless 
gauge boson associated with this U(l), and ^ and e* 
(i = 1, 2, • • ■ , I\Tj) represent the massless (4 component) 
fermions and their associated charges under U(l)^. Note 
that some of the ej may be zero, thereby allowing us to 
dial the number of free-streaming species. We consider a 
complex scalar particle, denoted by y, as the DM candi¬ 
date. 

Since our scenario assumes hidden sector dark matter, 
y is uncharged under the SM gauge interactions. The 
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only possible renormalizable interaction of \ with SM 
particles is of the Higgs portal form |x| 2 |H| 2 , where H is 
the complex scalar doublet that includes the SM Higgs 
particle h. This operator is expected to be present, and 
serves to ensure that the SM and hidden sectors are ini¬ 
tially in thermal and chemical equilibrium. However, our 
interest is in the scenario in which this coupling is not 
large enough to govern the annihilation of dark mat¬ 
ter that determines its relic abundance. Instead, the 
relic abundance is controlled by an additional interac¬ 
tion through which dark matter annihilates into DR. For 
that purpose, we introduce a massive spin-1 boson Z that 
couples to both x and ijji, which carry charges q x and (ji 
under the associated broken U(l)^. In our analysis, we 
will focus on the parameter range where x is lighter than 
the Z gauge boson by a factor of a few. 

The Lagrangian of our benchmark model therefore 
reads as 

C=C SM + (D M x)*(D /i x) - m 2 |x| 2 - *\x\ 2 \H\ 2 - l|x| 4 * 
+ hMi - \z^Z^ + l -m\Z 2 - 

(41) 

In this expression the ellipses represents the sector re¬ 
sponsible for breaking U(l)^ and generating the Z mass. 
We will not specify this sector explicitly as it is not rel¬ 
evant for our discussions. We only mention that all the 
particles in that sector are assumed to be significantly 
heavier than m x and m^, so that they do not affect the 
dynamics we are considering^] Note that, as a conse¬ 
quence of electroweak symmetry breaking, the parameter 
rh x in the Lagrangian is not equal to the x mass, m x , but 
is related to it as 

m x = ™x + f “ew , (42) 

where v EW = 246 GeV. The masses of the i/j are all set 
to zero, and are protected against quantum corrections 
by chiral symmetry. Finally, to ensure the stability of the 
dark matter particle x> we impose an exact 7Li symmetry 
under which x is odd and all the other fields are evenP] 

In this benchmark model, there is only one renormal¬ 
izable interaction, k|x| 2 |I?| 2 , that connects the dark sec¬ 
tor to the SM. The presence of such an interaction is 
expected to be a general feature of any theory in which 
the dark matter particle is a scalar, as there is no sym¬ 
metry that forbids it. Provided n > 10~ 6 , this coupling 


ensures that the SM and hidden sector are initially in 
thermal equilibrium at temperatures of order the weak 
scale. The requirement that annihilation to SM states 
does not play a significant role in setting the relic abun¬ 
dance constrains ac < 10 -2 . We will therefore focus on 
the regime 10 -2 > k > 10 -6 . We also choose to stay 
away from the region of parameter space such that the 
DM mass is close to half the Higgs mass, where annihi¬ 
lation to the SM is resonantly enhanced. 

The presence of interactions between the hidden sector 
and the SM leads to constraints on the theory and poten¬ 
tial signals. Upon electroweak symmetry breaking, the 
|x| 2 |H| 2 term leads to a 3-point interaction of the form 
/i|x| 2 - Provided m x < mh/2, this will generate a con¬ 
tribution to the invisible width of the Higgs boson given 

by 


r 


inv 

h 


m h ac 2 v 2 w „ 
16tt m 2 h Px ’ 



(43) 


The current experimental limit on the invisible branch¬ 
ing ratio from the LHC 7 and 8 TeV datasets stands at 
about 30% of the total Higgs width, = 4.07 x 10 -3 
GeV, from the vector boson fusion channel [651 EE] • The 
limit from the associated production channel is signifi¬ 
cantly weaker E3EE]. The bound on the invisible width 
is expected to improve to ~ 10% of the total width after 
300 fb" 1 at the 14 TeV LHC [H ED]. For m h > m x , 
these limits can be translated into the constraints on the 
parameter ac. At present we have ac < 2 x 10 -2 from 
the LHC 7 and 8 TeV runs. This limit is expected to 
improve to ac < 7 x 10~~ 3 after 300 fb -1 at the 14 TeV 
LHC. Stronger bounds can be obtained at future lep¬ 
ton colliders, such as the ILC or TLEP. According to the 
studies |7lH73j these machines can constrain the Higgs in¬ 
visible branching ratio to 0.2% ~ 1%, which corresponds 
to ac < 9.5 x 10~ 4 ~ 2.5 x 10 -3 . For m.h < 2 m x , on-shell 
Higgs decays into DM particles are kinematically forbid¬ 
den, and the hidden sector must be accessed through an 
off-shell Higgs, or through loop effects. The current col¬ 
lider limits are then very weak in the regime of interest, 
At < 10 — 2 , and are not expected to improve significantly 
even at future colliders [72] . While the h |x| 2 interac¬ 
tion can also give rise to signals in direct/indirect DM 
detection experiments, current experiments are not yet 
sensitive in the region At < 1(V 2 . 


B. Relic Abundance of DM 


2 We assume that the Lagrangian in Eq. © does not include 

kinetic mixings between these new U(l) gauge fields and those of 
the SM, such as , where B M is the SM hypercharge gauge 

boson. These can be forbidden by an internal charge conjugation 
symmetry carried by the hidden sector fields. 

3 Alternatively, such Z2 symmetry can accidentally emerge if the 

value of q x is such that it is not equal to any linear combination 

of the other U(l)^ charges with rational coefficients. 


Let us now analyze this benchmark model. The first 
step is to require that the relic abundance of x agree with 
the observed amount of cold dark matter. We focus on 
the region of parameter space where the DM mass lies 
between 5 GeV and 100 GeV, and where 2 m x < m 
Then the abundance of DM is governed by the annihila¬ 
tion process x + X* ~^ ^ ^ y ia the exchange of a virtual 
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Z in the s-channel, after the DM particles become non- 
relativistic. To lowest order in the nonrelativistic limit, 
the cross section for this process is given by 


where E^ represents the energy of the scattering ipi in 
the CM frame. Performing a thermal averaging, we ob¬ 
tain 




ffeff 

487n7i 2 


1 - 


4m* 


u x > 


(44) 


i =1 


4 • 45C(5) gj s f 5 
47t 2 27t 


= r col , (49) 


where 


N * 

~4 _ -4 -2 -2 

geff = 

i=l 


(45) 


and v x represents the speed of each annihilating \ i n the 
center-of-momentum (CM) frame. The presence of the 
v x suppression is an indication that the annihilation pro¬ 
cess proceeds through the p -wave channel. This can be 
understood from the fact that the initial state consists of 
two scalars and the intermediate state is spin-1. Letting v 
be the relative velocity of the annihilating x’s (v = 2v x ), 
and performing a thermal averaging, we obtain 


xx 


xifiip v ^ 


geff f 1 m ! \ 2 3T 

487T7T7. 2 I 4TO 2 ) m x 


(46) 


In this expression T is the temperature of the dark sec¬ 
tor and we have used the nonrelativistic relation ( v 2 ) = 
6 T/m x . We emphasize again that T is in general no 
longer equal to the temperature T of the SM gas, be¬ 
cause we are allowing for the possibility that the dark 
sector decouples from the SM well before before the \ 
particles freeze out. 

The relic abundance of x can be obtained from the 
Boltzmann equation, 

^ + 3 Hn x = ~(<r xx ,^ v) (n 2 - (n^) 2 ) . (47) 

This can be solved using the methods discussed in sec¬ 
tion [nil 


C. Kinetic Decoupling between DM and DR 

After freeze-out, the system of DM and DR continues 
to behave as a tightly coupled fluid until kinetic decou¬ 
pling occurs. This is in analogy to the photon-baryon 
plasma before recombination, and the DM-SM plasma in 
the context of the conventional WIMP scenario. Even 
though DM and DR have chemically decoupled, they are 
kept in kinetic equilibrium by the elastic scattering of \ 
with the fermions ipi, which proceeds through an off-shell 
Z in the t-channel. To leading order in the nonrelativis¬ 
tic limit, the cross section for \ + V’i —► X + is given 

by 


a 




M '2 ~2 E 2 7 
y 'i’i 

2n m\ 


(48) 


where the overall factor of 4 in front accounts for the two 
polarizations of ip, as well as the contribution from scat¬ 
tering with an anti-^j. The factor of 45<C(5)/47r 2 arises 
from thermal averaging. Note that T co i represents the 
rate for a x to experience a single collision with any one of 
the ipi. This rate is not, however, the same as the rate for 
the x plasma to thermalize with the ip bath. This is be¬ 
cause a single collision of a x particle with a ipi typically 
involves a momentum transfer of only (D(T), which is not 
large enough to significantly deflect the direction of the 

X , which carries a momentum 0(\Jm x T). Viewing the 
effect of N such collisions as a random walk with N steps, 
the x momentum typically changes by 0(VNT) after N 

collisions. Requiring that this change be 0(\Jm x T), we 
obtain that N ~ m x /T. It follows that the rate for the 
X plasma to thermalize with the ip bath is given by 

T coI 45C(5)5 e 4 ff T 6 

v 2o ™4 ™ 

l\ 71 ' 27T TO ~TO v 

z A 


Then, the temperature Tb below which the x plasma 
thermally decouples from the bath may be estimated as 


Tcoi 

IsT 


T=T d 


H 


T=T d 


(51) 


For the ranges of parameters considered in our bench¬ 
mark models, a quick estimate tells us that Tb ~ ^(l) - 
0(10) MeV. 

Above the kinetic decoupling temperature Tb, DM is 
in equilibrium with DR, and they form a tightly coupled 
fluid. Acoustic oscillations within this fluid have the ef¬ 
fect of erasing density perturbations on small scales. As 
a result, the temperature Tb determines the cutoff of the 
power spectrum, and sets a lower bound on the masses 
of the smallest halos m, 


m -= 105 (kSv)" S m ®' <52) 

In principle, this offers a separate way to probe these the¬ 
ories, independent of the CMB. At present, kinetic de¬ 
coupling temperatures up to about 10 keV can be probed. 
Unfortunately, in our scenario, Tb is too high for the cut¬ 
off of short-distance structures to be observable in current 
experiments. 

One possible generalization of our scenario that would 
lead to observable effects in the DM power spectrum in¬ 
volves allowing the mass of the mediator to lie below 














11 


the weak scale. It is easy to see that the WIMP mira¬ 
cle prediction for DM relic abundance continues to ap¬ 
ply, provided the DM candidate itself continues to have 
weak scale mass and couplings of 0(1) to the mediator. 
In order to have observable effects, a mediator of mass 
rhz < 0(100) MeV would be required. If the mediator 
mass is further pushed down to O(MeV), exchange of the 
light Z would give rise to sizable DM self-interactions, 
which can lead to observable effects in the DM power 
spectrum m • The resulting DM scattering cross section 
is velocity dependent, and can therefore be large enough 
to resolve small scale structure anomalies while remain¬ 
ing consistent with cosmological bounds, along the lines 
suggested in mmm- We leave a detailed study of 
this possibility for future work. 


D. Kinetic Decoupling of DM from the SM bath 

As we noted earlier, the size of the CMB signals de¬ 
pends on the temperature T = T = Tkd at which the dark 
matter sector kinetically decouples from the SM bath. 
For m x in the range between 10 GeV and 100 GeV, the 
leading processes maintaining the kinetic equilibrium are 
the elastic scattering of x off the W, Z , and SM fermions. 
These processes decouple at different temperatures, and 
Tkd is set by the last process that decouples. 

It turns out that at temperatures above ~ 10 GeV the 
dominant process is the scattering of DM particles off W 
bosons. At these temperatures, even if the W bosons 
are nonrelativistic and their number density begins to be 
Boltzmann suppressed, it is still more dominant than the 
competing scattering of x off a relativistic b quark, as the 
latter suffers from a small Yukawa coupling. Let’s first 
consider the scattering of a relativistic x off a nonrela¬ 
tivistic W. The tree-level cross section in this limit is 
given by 


2 2 
_ k m z w 

4?r m A h ' 


(53) 


Equating the scattering rate = n w (a w v) with the 
Hubble expansion rate H gives an estimate for the tem¬ 
perature at which this process decouples, where n w rep¬ 
resents the equilibrium number density of W bosons, and 
the relative speed v between the scattering \ and W is 
unity in the limit under consideration. We thus obtain 


Tkd 


m w 



9\/l0 Mp rn^y n 2 

167T 4 m h m\ v / g * 


-l 


= rn w 


8.5 + log 


k 2 10 V 

io - 10 y/g*)_ 


(54) 


In obtaining this expression, we have assumed that g * 
g*, so that the energy density is dominated by the SM 
degrees of freedom. For a similar process with the Z bo¬ 
son rather than the W, every appearance of m w above 


should be replaced by m z , and the argument of the log¬ 
arithm is reduced by a factor of 2 because of the fewer 
degrees of freedom associated with the Z. Therefore, the 
corresponding decoupling temperature is proportional to 
m z with the same proportionality factor (neglecting the 
change in the argument of the logarithm), so the Z de¬ 
couples at a slightly higher temperature than the W. 

On the other hand, in the limit where the x and W are 
both nonrelativistic, the scattering cross section is given 
at tree level by 


u w 


2 4 

47 t ( to x + m w ) 2 ml 


(55) 


To calculate the rate n w (a w v), we need to know (v). 
For a single nonrelativistic particle of mass m obeying the 
Maxwell-Boltzmann distribution with temperature T, an 
elementary integral gives (v) = yj8T/(nm ). We have a 
two-body system instead, but the relative speed v in the 
lab frame is the same as that in the CM frame, where 
the two-body problem reduces to a one-body problem 
with the reduced mass p xw = m x m w / (m x + m w ). We 
thus have (v) = \J 8T/(7r p xw ). The kinetic decoupling 
temperature then becomes 


Tkd 


-l 


8.9 + log 


A*rw 


10 





(56) 

At temperatures below ~ 10 GeV, scattering of DM off 
SM fermions becomes increasingly important, and even¬ 
tually dominates. At these temperatures the DM parti¬ 
cle x is nonrelativistic, while the SM fermions may be 
relativistic or nonrelativistic. In the case where the scat¬ 
tering SM fermion / is relativistic, the tree-level cross 
section is given by 


_ k 2 m 2 p 2 
8nm x ml 


(57) 


where p is the magnitude of the 3-nromenta of the scatter¬ 
ing particles in the CM frame, and higher order terms in 
p/m x have been neglected. In this leading nonrelativis¬ 
tic approximation, the thermal average (p 2 ) in the CM 
frame is the same as that in the lab frame, as the differ¬ 
ence between the two frames itself is an 0(p/m x ) effect. 
Therefore, we evaluate (p 2 ) and the number density n/ 
simply in the lab frame in the standard way, obtaining 


(P 2 ) = 


15C(5) 2 
C(3) 


»/ = «c^r s 


(58) 


where N c is the number of colors carried by the fermion 
species /. Notice, however, that the relevant reaction 
rate to maintain a kinetic equilibrium is not simply given 
by rif(ijfv), because a single collision of a heavy parti¬ 
cle x with a relativistic particle / with p m x hardly 
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SM-DS Kinetic Decoupling Temperature 



Note that the “reference value” of k used here is 10 -3 , 
which differs from Eqs. (541, (56), and (59), where the 
value used was 10 -5 . This reflects the fact that this pro¬ 
cess is relevant only if k is sufficiently large, even for the 
b quark. 

In Fig. [4j we plot the kinetic decoupling temperature 
against the dark matter mass, for several different values 
of the coupling n. We see that for the range of values of 
k that can be probed in current and future collider ex¬ 
periments, the kinetic decoupling temperature lies below 
several GeV and above several hundred MeV. In particu¬ 
lar, the kinetic decoupling occurs before a drastic change 
in the number of relativistic SM degrees of freedom due 
to QCD phase transition. 


FIG. 4. The kinetic decoupling temperature Tkd between the 
SM and hidden sectors as a function of the DM mass todm 
for three values of k: klhci4 = 7 x 10 -3 , kilco2 = 9.5 x ICG 4 
and k = 5 x 10 -6 . 


changes the direction of the x momentum. Since the \ 
momentum is on average of order ^/m x T, and a typ¬ 
ical momentum transfer via a single collision is of or¬ 
der p ~ T, the number of collisions N required to al¬ 
ter the x momentum by an 0(1) fraction is given by 
y/NT ~ yjm x T , i.e., N ~ m x /T. The relevant reaction 
rate to maintain a kinetic equilibrium, therefore, is given 
by Tf = nf(crfv)/N. Then, from T/ ~ H, we obtain an 
estimate for the kinetic decoupling temperature: 


/ 167T 9 / 2 1 to 3 y/g„ V 

\ y 135v / 5C(5) N c m 2 f M P k 2 J 


~ 4.3 GeV x 


( 3 m 3 m 2 10- lo V^Y 
yN c (10 GeV ) 3 to 2 k 2 10 j ' 

(59) 


where mb is the 6 -quark mass. 

If the SM fermion is also nonrelativistic, the tree-level 
cross section in the leading nonrelativistic approximation 
is given by 




An(m x + to /) 2 


(60) 


For thermal averaging, we use (v) = yj8T / (tt p x f ), sim¬ 
ilarly to what we did below Eq. (551. Furthermore, the 


number of collisions N required to randomize the motion 
of x is given by '/NyJmfT ~ yjm x T , i.e., N ~ to x /to/. 
Then, from rif(<jfv) /N ~ H , we obtain 


T kd ~ TO/ 


= TO/ 


log ( ' V5 • N t: ■ Mp 

\^27T 9 / 2 m h m\m\ \ to/ / yfg* J 


5.2 + log |V< 



E. Signals 


Scattering (\j/ + A^) 


- [1,14-7 X 10 3 

- ailco 2=9.5 x 10 

- K=5 x 10" 6 


■x; 


PIanck2015(TT)+LowP+BAO 


15 20 30 

m dm (GeV) 


a= 

Z 

<1 0.19 

0.15 

0.12 

0.10 


Free Streaming (if/) 


50 70 100 


- *LHC14=7 X 10 

- *ILC02=9.5 x 10" 

- K =5 x 10 -6 


X 


Planck2015(TT)+LowP+BAO 


-0.87 

0.69 

-0.55 

0.46 


10 15 20 30 

m dm (GeV) 


50 70 100 


FIG. 5. AlVefi and At as a function of todm- It is as¬ 
sumed that > Todm- Three values of k are considered: 
klhci4, «ilco2 and k = 5 x 10~ 6 . Also shown are the 2015 
Planck results: the central value (Green dashed line) and the 
2 a constraint (Orange dashed line). 


Given the dark matter mass and the value of the cou¬ 
pling constant k. we can determine the temperature T kd 
at which the dark sector kinetically decouples from the 
SM. We can then use Eq. ([3]) to obtain AAT e ff and 
Eq. (12) to obtain the angular shifts of the CMB peaks. 
In our benchmark model, g* and < 7 * k d are given by 


9* ~ 2 + 2 N $ ■ 


5 * kd =n+ -N^ . 


(62) 
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We assume that the Z boson is already nonrelativistic 
when the SM and hidden sectors decouple. Then n = 2 
if the DM candidate is already nonrelativistic when the 
SM and dark sectors decouple, and n = 4 otherwise. 

The results for AO = 1 are shown in the upper panel 
of Fig. [5] for three different values of k. The smallest 
of the values studied, k = 5 x 10 —6 , corresponds to the 
minimum value of k that can ensure that the SM and 
hidden sector are in thermal equilibrium at or above the 
weak scale. The other two values studied correspond to 
the limits that can be placed on n at the LHC and at a 
lepton collider. In the figure we have plotted the current 
limits on A7V e ff from the Planck experiment. Although in 
the region of light DM masses that may be probed at the 
14 TeV LHC the predicted A N e g are above the central 
value from Planck data fit, they are still well within 95% 
C.L., and are large enough to be detected in upcoming ex¬ 
periments. It is also important to note that these bounds 
are not immediately applicable to this model since in the 
standard Planck analysis, it is assumed that the contri¬ 
bution to A N e g is free streaming. A fresh analysis that 
relaxes this assumption by including scattering radiation 
is required to determine the current limits on this sce¬ 
nario. 

We see from Fig. [5] that A N e g > 0.2 in this class of 
models for the entire range of DM masses. Then, based 
on the discussion in Section II we expect that upcoming 
experiments will offer the possibility of distinguishing this 
scenario from one where A N e g is the same, but the DR 
free streams. For the purposes of comparison, in the 
lower panel of Fig. [5] we have plotted A N e g and A£ for 
the same model, but without the U(l) gauge symmetry. 
The DR now consists of only one massless fermion species 
Since the gauge boson A M is now absent, the DR free 
streams rather than scatters. We see from the figure that 
although the magnitudes of AJV e g and A£ are comparable 
in size to the benchmark model, the sign of A£ is opposite 
in sign. It is this difference that we expect will help 
distinguish between free streaming and scattering DR. 


V. CONCLUSIONS 

In this paper, we have explored a scenario where the 
DM candidate is part of a hidden sector whose parti¬ 
cles carry no charges under the SM gauge groups, and 
whose couplings to the SM states, though present, are 
small. The abundance of DM is assumed to be deter¬ 
mined primarily by annihilation to massless states that 
lie within the hidden sector. Then, if we further assume 
that the weak scale is the only mass scale in the hid¬ 
den sector, so that the mass of the DM particle and its 
annihilation cross section are both set by this scale, the 
observed abundance of DM can be naturally explained. 
This framework is motivated by an alternative realiza¬ 
tion of thermal WIMP DM paradigm that is naturally 
compatible with limits from direct, indirect and collider 
searches for DM thus far. 


A robust consequence of this framework is the exis¬ 
tence of DR, associated with the massless states in the 
hidden sector. The contribution of this DR to the energy 
density of the universe during the era of recombination 
epoch manifests itself observationally as a contribution to 
the effective number of neutrino species, AA e ff. In addi¬ 
tion, massless particles constituting DR may or may not 
interact with one another. We determined the shift in 
the locations of the CMB peaks, A£, as a function of the 
number of free streaming DR species, A N^ e , and scat¬ 
tering species, A7V|g. att . We found that free streaming 
and scattering species shift the peaks in opposite direc¬ 
tions, so that by combining this effect with the measure¬ 
ment of the total A N e s mentioned above, we can sepa¬ 
rately determine A N^' 6 and A N^ att . We also calculated 
the corrections to the amplitudes of the scalar and ten¬ 
sor modes of the CMB arising from the presence of DR, 
and showed that the sign of the correction depends on 
whether the DR scatters or free streams. 

We have found that provided the hidden sector is ini¬ 
tially in thermal equilibrium with the SM degrees of free¬ 
dom at temperatures at or above the weak scale, there is 
a robust prediction for a lower bound on A N e g of about 
0.02. This is large enough to be observed by future CMB 
Stage-IV experiments. In the scenario where the DR 
has self-interactions, assuming naturalness, A N e g is ex¬ 
pected to lie well above this lower bound, making it large 
enough to be observed in upcoming experiments such as 
CMBPol. We further constructed a simple model that 
realizes this scenario. In our model, the SM and hidden 
sector are initially kept in thermal equilibrium through 
the Higgs portal. We have determined the size of the 
AiV e ff signal, and shown that it is large enough to be 
detected in upcoming experiments. These experiments 
are also expected to be sensitive enough to the A£ sig¬ 
nal to allow the possibility of distinguishing between free 
streaming DR and scattering DR. In addition, we find 
that there are regions of parameter space where invisible 
decays of the Higgs into hidden sector states can be used 
to probe this model at colliders. 

An important point to reiterate is that the bounds on 
AA e ff derived from the standard analyses, such as by 
Planck, may not be directly applicable to our scenario 
where the DR may be self-scattering. A dedicated anal¬ 
ysis parametrized by AN^ e and AfVJg att instead of a 
single AN e s is required to determine the exact limits on 
this scenario. 

The general possibility that dark matter may arise 
from a sector hidden from the visible matter of the SM 
has drawn considerable interest recently, due in no small 
part to the increasing experimental limits on the conven¬ 
tional WIMP DM paradigm. Our work highlights that 
in this scenario which is challenging or perhaps even im¬ 
possible for conventional DM searches, CMB experiments 
may be able to shed light on the nature of the DM sector, 
which is worthy of further exploration. 
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